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ABSTRACT 


The  mechanics  of  separation  of  a  thin  interfacial  liquid  layer  trapped  between  two 
parallel  surfaces  was  studied  in  a  controlled  manner.  Different  liquid  viscosities,  layer 
thicknesses  and  separation  velocities  were  used  to  investigate  the  extensional  behavior 
and  determine  its  dependence  on  viscous  fingering,  and  Capillary  number.  Force, 
displacement  and  time  data  have  been  recorded  for  all  experimental  runs.  Qualitative 
visual  data  have  also  been  recorded  to  corroborate  the  trends  in  the  onset  of  viscous 
fingering  based  on  a  simple  interfacial  stability  analysis.  The  quantitative  data  has  been 
used  to  generate  force-displacement  plots  of  the  separation.  The  results  of  this  work 
provide  useful  fundamental  insight  into  the  mechanics  of  this  novel  problem. 
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I.  INTRODUCTION 


The  original  studies  on  the  instability  of  the  interface  between  two  liquids  were 
carried  out  in  the  late  1950s,  when  water  was  pumped  down  into  an  apparently  ‘dry’  oil 
field  to  push  up  the  remaining  oil.  The  belief  was  that  the  oil  being  lighter  than  water 
would  rise,  although  in  reality  the  interface  between  the  water  and  oil  was  unstable  and 
fingers  of  water  penetrated  up  through  the  ground  leaving  the  oil  down  below.  It  was 
shown  by  Taylor  (1958),  and  Saf&nan  and  Taylor  (1958),  in  a  now  classic  study  that  the 
phenomenon  could  be  explained  by  a  simple  analysis  based  on  Darcy’s  law  for  flow  of  a 
viscous  fluid  through  a  homogeneous  porous  medium. 

The  study  of  the  motion  of  fingers  of  a  less  viscous  fluid  through  a  more  viscous 
fluid  in  a  Hele-Shaw  cell  (i.e.  in  a  narrow  gap  between  parallel  plates.)  is  a  paradigm  in 
applied  mathematics.  It  has  raised  questions  at  all  levels.  By  choosing  the  proper 
equations  of  motion,  it  can  be  solved  analytically  or  numerically.  Considerable 
experimental  evidence  has  also  been  gathered  for  various  scenarios.  Continuing 
disagreement  with  experiments  leads  to  modifications  of  the  assumptions  and 
approximations  employed  by  the  theory  until  good  agreement  is  reached  to  help 
understand  the  phenomenon  (Saffman,  1991). 

It  has  been  shown  that  for  a  channel  of  given  aspect  ratio  (width/thickness  of  the 
layer),  it  was  foimd  for  the  case  when  the  viscosity  of  the  displacing  fluid  is  negligible 
that  the  shape  and  width  of  the  fingers  depend  only  on  the  Capillary  number  (Ca).  A  way 
to  measure  the  ratio  of  viscous  forces  to  surface  tension  is  through  the  use  of  the 
Capillary  number  (Ca  =  pU/a).  If  the  capillary  number  remains  small  (10~^-10’^),  the 
surface  tension  effects  are  expected  to  dominate  the  flow.  The  combination  of  the  liquid 
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viscosity  and  interface  velocity  could  be  chosen  in  such  a  way  that  the  capillary  number 
would  be  in  the  above  limits.  Also,  by  keeping  the  aspect  ratio  large  (D/lo  »1),  the  fluid 
mechanics  in  the  layer  could  be  restricted  to  essentially  one-dimensional  behavior. 

In  this  experiment,  the  approach  to  the  problem  is  different.  A  tensile  force  is 
applied  to  a  layer  of  liquid  trapped  between  two  circular  parallel  disks  and  the  force- 
displacement  behavior  is  studied  vis-a-vis  the  onset  of  viscous  fingering.  A  -wide  range 
of  viscosities,  separation  velocities  and  initial  liquid  layer  thicknesses  are  explored. 

A  variation  of  the  Saffinan-Taylor  instability  analysis  will  be  used  to  determine 
the  neutral  stability  curve  and  the  most  dangerous  wavelength  that  grows  into  fingers,  as 
well  as  the  number  of  corrugations  (fingers)  that  form.  Through  experimental  analysis  of 
the  data,  the  interface  velocity  and  the  Capillary  number  at  the  point  of  viscous  fingering 
can  be  calculated  with  extreme  accuracy  and  repetition. 

Quantitative  force-displacement  data  and  qualitative  visual  data  have  been 
gathered  for  this  experimental  study.  Both  forms  of  data  have  been  used  to  make  some 
useful  deductions  of  the  mechanics  of  separation  of  thin  interfacial  liquid  layers. 
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II.  BACKGROUND 


A.  SURFACE  TENSION 

1.  Theoretical  Background 

The  effects  of  surface  and  interfacial  tensions  give  rise  to  so  many  phenomena 
observed  in  liquid  behavior  that  the  complex  physical-chemical  interactions  involved  are 
taken  for  granted  and  not  all  of  which  are  understood  even  today.  The  liquid  state  itself  is 
composed  of  molecules  in  motion  that  are  kept  relatively  close  to  each  other  by  attractive 
Van  der  Waals  forces.  However,  a  principal  method  of  analysis  of  problems  of 
interfacial  effects  rests  upon  the  assumption  that  a  mean  molecular  field  can  describe  the 
liquid.  It  is  assumed  possible  to  define  an  element  of  the  liquid  that  is  small  compared  to 
the  range  of  intermolecular  forces  but  large  enough  to  contain  a  sufficient  number  of 
molecules.  This  approximation  implies  that  on  average  the  attractive  force  on  any 
molecule  in  the  liquid  is  the  same  in  all  directions  giving  to  the  liquid  its  fluid 
characteristics  (Probstein,  1989). 

2.  Liquid-Gas  Interface 

At  a  liquid-gas  interface  although  the  molecules  are  fi-ee  to  move  in  the  liquid, 
their  motion  is  far  more  restricted  than  in  a  gas  where  there  is  a  little  attraction  between 
the  molecules.  Therefore  the  attraction  between  the  liquid  molecules  will  prevail  and 
prevent  the  liquid  molecules  fi-om  escaping  into  the  gas.  As  a  result,  the  liquid  molecules 
at  the  surface  are  attracted  inward  and  normal  to  the  liquid-gas  interface,  which  us 
equivalent  to  the  tendency  of  the  surface  to  contract. 

If  the  interface  is  curved,  a  mechanical  balance  shows  that  there  is  a  pressure 
difference  across  the  interface.  The  pressure  is  higher  on  the  concave  side.  Two  surface 
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tension  forces  balance  the  pressure  increase  in  the  interior.  For  an  arbitrary  curved 
interface  whose  principal  radii  of  curvature  are  Ri  and  R2  and  have  a  surface  tension  ct. 
A  force  balance  normal  to  the  surface  will  show  that  the  pressure  differential  is 

A  f  ^  0 

The  previous  equation  is  known  as  the  Young-Laplace  Equation.  For  a  special  case  of  a 
spherical  bubble  or  drop  in  immiscible  liquids,  a  well-known  result  is  obtained,  where  r  is 
the  bubble  or  drop  radius. 


Lp  = 


2(7 


(2) 


When  a  drop  is  placed  on  a  plane  solid  surface,  it  will  be  in  contact  not  only  with 
the  surface  but  also  vidth  a  gas  such  as  air.  The  liquid  may  spread  freely  over  the  surface, 
or  it  may  remain  as  a  drop  with  a  specific  angle  of  contact  with  the  solid  surface.  A  force 
balance  would  involve  both  surface  tension  (cr)  and  contact  angle  (0).  If  the  contact  angle 

is  less  than  90°,  the  liquid  is  said  to  wet  the  solid;  if  the  contact  angle  is  greater  than  90°, 
the  liquid  is  nonwetting  (White,  1994). 
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3.  Capillaiy  Number  (Ca) 

A  dimensionless  quantity  which  plays  a  significant  role  in  determining  the  effect 
at  the  interface  of  the  liquid-gas  boundary  is  the  capillary  number  (Ca).  The  capillary 
number  measures  the  ratio  of  the  viscous  force  to  surface  tension  force  and  is  defined  by 
A  small  capillary  number  indicates  that  the  viscous  forces  are  not  dominating  the  flow, 

_  viscous  force  uU 

Ca  =  — — - ; - =  .c —  (3) 

surface  tension  force  a 

the  surface  tension  acts  as  a  stabilizing  factor  in  examining  the  effects  of  viscous 
fingering. 

B.  VISCOUS  FINGERING 

1.  Theoretical  Background 

An  important  instability  phenomenon  that  has  a  widespread  importance  in  this 
experiment  is  viscous  fingering.  This  is  essentially  an  instability  of  bounded  creeping 
flows,  such  as  appear  in  a  porous  medium  or  a  Hele-Shaw  cell,  when  one  fluid  is 
displaced  by  a  second  fluid  of  lower  viscosity.  Differential  pressures  imposed  at  the  open 
boundaries  could  be  the  driving  mechanism  for  the  flow.  The  fundamentals  of  creeping 
flows  and  Hele-Shaw  flows  are  first  briefly  reviewed  before  applying  them  to  the  viscous 
fingering  relevant  to  this  study. 

2.  Creeping  Flows 

Some  approximate  solutions  of  the  Navier-Stokes  equations  for  a  limited  case 
when  the  viscous  forces  are  considerably  greater  than  the  inertia  forces  can  be  made. 
Since  the  inertia  forces  are  proportional  to  the  square  of  the  velocity  whereas  the  viscous 
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forces  are  only  proportional  to  its  first  power,  it  is  easy  to  appreciate  mathematically  a 
flow  for  which  viscous  forces  are  dominant.  Under  these  conditions  the  Reynolds 
number  is  very  small  and  the  inertia  terms  can  be  simply  omitted  from  the  equations  of 

motion  as  a  first  order  approximation.  It  can  be  seen  that  the  incompressible  form  of  the 
Navier-Stokes  equations  given  by 

Dsv 

-gradp  +  //V  %  (4) 

with  the  inertia  terms  neglected  become 

gradp= 

The  same  boundary  conditions  as  the  full  Navier-Stokes  equations  must  be  introduced  to 
the  system  of  equations,  namely  those  expressing  the  absence  of  fluid  slip  at  the  walls. 

3.  Hele-Shaw  Flows 

A  solution  of  the  three-dimensional  equations  of  creeping  motion  can  be  obtained 
for  the  case  of  flow  between  two  parallel  flat  walls  separated  by  a  small  distance  2h.  If  a 
cylindrical  body  of  arbitrary  cross-section  is  inserted  between  the  two  plates  at  right 
angles  so  that  it  completely  fills  the  space  between  them,  the  resulting  streamlines  is 
identical  with  that  in  potential  flows  about  the  same  shape.  It  can  be  shown  that  the 
solution  for  creeping  motion  then  possesses  the  same  streamlines  as  the  corresponding 
potential  flow  (Schlichting,  1979). 
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We  select  a  system  of  coordinates  with  its  origin  in  the  center  between  the  two 
plates,  and  make  the  x,  y-plane  parallel  to  the  plates,  the  z-axis  being  perpendicular  to 

them  as  in  Figure  2.  The  body  is  assumed  to  be  placed  in  a  stream  of  velocity  U,  parallel 
to  the  x-axis. 

From  the  continuity  equation,  assuming  two-dimensional,  incompressible,  fully 
developed  flow 


f-" 


(6) 


In  addition,  the  x-momentum  equation  simplifies  to 

dx  ^  dz^ 


(7) 


Integrating  with  respect  to  z  and  applying  the  boundary  conditions,  5u/5z  =  0  @ 
z  =  0  and  u  =  0  @  z  =  ±h/2  leads  to  a  velocity  profile 


fi  dx  2\  4  J 


(8) 
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From  the  velocity  profile  the  average  velocity  is  obtained 


12//  dx 

4.  Viscous  Fingering 

The  Saf&nan-Taylor  instability  at  an  interface  could  arise  when  two  fluids  of 
different  viscosity  are  pushed  by  a  pressure  gradient  through  a  Hele-Shaw  cell.  The 
instability  is  in  many  ways  similar  to  the  Rayleigh-Taylor  instability;  the  primary 
difference  is  that  the  equilibrium  state  in  the  former  is  a  dynamic  one,  in  which  the 
interface  between  the  two  fluids  is  moving  rather  than  stationary  (Faber,  1997). 
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Figure  3.  Side  View  of  Radial  Hele-Shaw  Cell. 


Here  we  closely  follow  and  adapt  the  development  by  Faber  (1997).  Suppose  the 
cell  is  horizontal,  consisting  of  two  parallel  radial  disks,  gravitational  effects  can  be 
ignored.  There  is  a  pressure  gradient,  which  drives  the  fluid  in  the  radial  (+r)  direction 
with  some  uniform  velocity  (Ur).  At  equilibrium,  the  interface  between  the  two  fluids  is 
a  concentric  circle  R  =  Uft.  Where  r  <  Urt,  the  fluid  viscosity  is  denoted  by  p';  where  r 
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>Ur  t,  the  viscosity  is  given  by  |t.  Using  the  Navier-Stokes  equations  in  cylindrical 
coordinates,  the  pressure  gradient  needed  to  maintain  this  motion  of  the  two  fluid  regions 
are  given  by 

4?'  llM'U,  4?  \2fiU^  (10) 

dr  dr 

Where  Ur  is  the  average  interface  velocity  and  d  is  the  disk  separation  distance.  The 
pressures  p'  and  p  are  not  necessarily  equal  at  the  interface,  because  the  interface  may  be 
curved  in  the  z  direction. 


+  P'  =  -^-^^^[r-Ud)  +  Po 


(11) 


Now  suppose  that  the  interface  is  perturbed  in  a  way  that  at  time  t  the  nominally  circular 
interface  lies  at  r  =  R,  where 


R  =  (12) 

k 


Interface 

Figure  4.  Plan  View  of  Radial  Hele-Shaw  Cell. 
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There  must  be  some  corresponding  perturbation  in  p'  and  p,  and  it  must  have  the 
same  periodicity  in  the  r  direction.  Any  perturbation  can  be  expressed  as  an  infinite  sum 
of  wavelike  functions.  Since  the  perturbation  cannot  effect  the  pressure  at  large  distances, 
these  pressure  terms  have  the  form 


P=-  -^72  '  {r -  Urt)  +  p,^Ce 


k{r-U,l) 


(13) 


p'  =  - 


Up'U, 


where  k  —  ItzIX,  and  Q'  and  C  are  constants  to  be  determined  by  boundary  conditions  at 
the  interface.  The  boundary  conditions,  applicable  in  each  case  at  r  =  R,  and  linearized  by 
omission  of  the  higher  order  terms  yield 


(14) 


12//'  dr  \lp  dr  'dt 


(15) 


which  corresponds  to 


d^k  d( 
12//'  12//  dt 


(16) 


The  pressure  differential  due  to  the  interfacial  surface  tension  is 


/9^  p 

W  ^  -1.2  nly  ^ikre 

P  P  ^  ^2  “  ^ 


(17) 
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The  corresponding  growth  rate  Sk  can  be  found  to  be 


1 


ak^d 


3  J2 


Ct  ^  //  +  //'  12 


-u,k{fi-A 


(18) 


It  may  be  noted  that  the  interface  is  stable  for  all  values  of  k  when  p'  <  p.  When 
p  <  p',  i.e.  a  viscous  fluid  is  being  displaced  by  a  less  viscous  one,  the  interface  can 
become  imstable.  The  neutral  stability  curve  corresponding  to  zero  growth  rate  shows 
that  it  is  marginally  stable  for  the  critical  wavenumber  k  =  kc  where 


UuXfi'-M) 

ad^ 


(19) 


So  for  wavelengths  oo  >  X  >  X,c  the  interface  is  imstable.  The  perturbation  that  grows  the 
fastest  is  when  Sk  is  a  maximum,  which  yields  the  most  dangerous  wavelength  as 


k  = 

max  ^ 


A  =  A  J3 

'‘'max  '‘'c  » 


or 


(20) 


(21) 


This  corresponds  to  a  wavelength  (A.max)  which  is  the  fastest  growing  disturbance  and 
hence  the  one  most  likely  to  develop  into  fingers 


'^max  ~ 


U{ju'  -  m) 


(22) 
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The  largest  value  of  the  wavelength,  which  is  consistent  with  the  circular  geometry  of  the 
flow  cell  of  radius,  r,  is 


'^/mgest  ~ 


(23) 


Then  for  a  given  combination  of  interface  velocity,  gap  separation,  surface  tension  and 
fluid  viscosities,  the  most  dangerous  wavelength  can  be  foxmd  and  the  number  of  fingers 
that  consequently  appear  can  be  calculated  from 


N  = 


Ikt 


(24) 


Figure  5  is  a  typical  plot  of  the  growth  rate  as  a  function  of  wavelength  for  a 
given  separation  velocity.  The  plot  clearly  shows  that  the  viscosity  of  the  fluid  affects 
growth  rate  and  its  magnitude  increases  with  increasing  viscosity.  It  can  also  be  said  that 
the  maximum  growth  rate  is  not  only  larger,  but  the  peak  is  also  more  pronounced,  for 
higher  viscosities. 
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C.  PREVIOUS  STUDIES 

The  paper  by  Paterson  (1981)  is  the  only  earlier  study  to  have  dealt  with  radially 
driven  flow  in  a  circular  Hele-Shaw  cell.  An  approximate  equation  for  the  growth  of 
viscous  fingers  was  proposed.  To  model  the  initial  growth  of  fingers,  Darcy’s  Law  for 
governing  the  velocity  of  flow  in  a  porous  medium  or  Hele-Shaw  was  used.  An  equation 
for  critical  wavelengtihi  was  obtained,  and  it  was  shown  that  if  the  circxunference  of  the 
injected  bubble  was  less  than  the  critical  wavelength  the  interface  was  stable  and  no 
viscous  fingering  was  present.  In  the  review  by  Saffinan  (1991)  the  selection  mechanism 
and  stability  of  fingers  and  bubbles  in  the  Hele-Shaw  Cells  were  discussed.  The  main 
focus  was  to  better  xmderstand  the  uncertainties  with  respect  to  anomalously  thin  fingers 
and  the  nature  of  the  finite  amplitude  of  the  instabilities.  He  suggested  that  the 
uncertainty  in  the  boundary  conditions  at  the  moving  interface  may  be  the  cause  for  the 
difference  in  the  theoretical  and  experimental  results.  It  was  also  noted  that  for  unsteady 
fingering  the  problem  is  not  yet  resolved  and  many  open  questions  remain  especially 
since  unsteady  flows  form  singularities  in  finite  time. 

In  a  recent  experimental  study  by  Maxworthy  (1989),  the  study  of  the  stability  of 
the  circular  interface  in  a  Hele-Shaw  when  air  displaces  a  viscous  oil  has  been  used  to 
estimate  the  critical  wavelength  for  large  capillary  numbers  (Ca).  A  number  of  questions 
have  arisen  about  determining  the  exact  boundary  conditions  to  be  applied  at  the 
interface,  which  was  slightly  modified  in  this  experimental  study.  The  experimental 
results  had  a  tendency  to  approach  a  wavelength  that  is  a  constant  multiple  of  the  gap 
width  for  values  of  Ca  greater  than  unity.  The  reasons  for  this  discrepancy  are  due  to  the 
non-zero  finite  nature  of  the  Reynolds  number,  and  three-dimensional  effects. 
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III.  OBJECTIVES  OF  THE  PRESENT  STUDY 


The  objective  of  the  present  study  is  to  obtain  some  fundamental  insight  into  the 
onset  of  viscous  fingering  subjected  to  tensile  forces.  Unlike  previous  studies  in  this 
field,  the  main  focus  is  the  mechanics  of  separation  of  a  thin  liquid  layer  trapped  between 
two  parallel  surfaces,  which  are  being  pulled  apart  at  a  constant  velocity  normal  to  the 
liquid  layer.  There  is  no  shear  or  compressive  forces  acting  on  the  sample.  A 
mathematical  equation  is  derived  assuming  a  Hele-Shaw  relationship  for  the  onset  of 
viscous  fingering,  known  as  the  critical  wavelength.  A  similar  relationship  will  also  be 
developed  for  the  maximum  growth  rate  of  the  viscous  fingers,  known  as  the  maximum 
wavelength,  which  will  then  determine  the  number  of  fingers  produced.  The  velocity  of 
the  interface  and  separation  of  the  disks  will  be  investigated.  Experiments  have  been 
conducted  and  will  be  correlated  to  the  mathematical  equation  to  the  onset  and  mayimum 
growth  of  the  viscous  fingers.  The  work  is  experimental  in  nature  and  is  aimed  to  gather 
both  qualitative  and  quantitative  data.  Experiments  have  been  planned  and  conducted  in  a 
way  that  some  fundamental  insights  about  the  effects  of  critical  wavelength  interface 
velocity,  separation  distance  and  fluid  viscosity.  This  study  also  aims  to  analyze  the 
different  competing  regimes  in  the  fluid  deformation  process. 
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IV.  EXPERIMENT 


A.  SELECTION  OF  THE  MODEL  SYSTEM 

Since  this  experimental  study  involves  a  liquid  layer  in  tension  (normal  to  the 
plane  of  the  layer),  the  experimental  apparatus  was  constructed  to  apply  a  purely  tensile 
force  in  the  vertical  direction.  The  initial  separation  of  the  plexi-glass  disks  in  which  the 
liquid  layer  resides  is  about  100  pm.  The  stage  that  contained  the  plexi-glass  disks  was 
constructed  to  ensure  that  the  mating  surfaces  of  the  disks  were  parallel  to  one  another. 
The  plexi-glass  disks  were  also  developed  to  ensure  parallelism  of  their  mating  own 
surfaces.  The  vertical  travel  of  the  positioning  actuator  was  aligned  to  prevent  any  shear 
force  to  develop  during  the  experimental  procedure.  The  alignment  of  the  vertical 
actuator  and  plexi-glass  disks  were  perfectly  perpendictilar  to  the  direction  of  the 
separation  movement.  Figure  6  shows  the  experimental  apparatus,  the  supporting 
structure  was  constructed  of  high  grade  aluminum  and  mounted  on  a  permanent  rigid 
support  platform  to  prevent  any  deviation  from  the  normal. 


Figure  6.  Experimental  Set-Up,  Apparatus. 
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The  experimental  set-up  was  designed  to  gather  both  numerical  and  visual  data. 
One  CCD  camera  was  placed  in  the  set-up  to  provide  this  required  data.  Since  one  of  the 
objectives  of  the  experiment  was  to  determine  onset  of  viscous  fingering,  the  CCD 
camera  was  positioned  to  record  the  movemehts  of  the  liquid  layer  from  the  bottom  of  the 
plexi-glass  disks.  This  camera  was  mounted  on  the  vertical  moving  stage  (to  move 
together  with  the  disks)  to  prevent  the  focal  length  from  changing  during  the  operation  of 
the  experiment  (Figure  7).  The  CCD  camera  was  connected  to  monitors  for  viewing  and 
recording. 
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The  computer  was  used  to  run  and  store  data  for  latter  use  and  analysis  (Figure  8). 
A  software  program  written  for  this  particulEir  experiment  controlled  the  vertical 
movement  of  the  actuator.  The  software  program  enabled  us  to  control  the  actuator  by  the 
computer  and  to  maximize  the  accuracy  of  the  experiment.  From  the  control  unit, 
displacement  of  the  discs,  time  and  from  the  load  cell,  the  force  information  was 
recorded. 

The  experimental  set-up  was  designed  to  give  numerical  data  and  visual  data;  the 
visual  data  could  be  viewed  and  recorded  via  the  CCD  camera  (Figure  9).  The  numerical 
data  was  achieved  through  the  use  of  a  load  cell,  which  was  mounted  on  top  of  the  upper 
plexi-glass  disk  (Figures  10  and  11).  The  load  cell  is  connected  to  an  amplifier  and 
multi-meter.  The  load  cell  was  used  to  measure  the  tensile  load  acting  of  the  fluid;  the 
data  was  then  stored  into  the  computer  for  further  analysis  and  plotting. 
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Figure  11.  Close-up  view  from  the  top. 
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B.  EXPERIMENTAL  SET-UP 


This  section  introduces  each  component  of  the  experimental  apparatus  and 
includes  a  brief  description  of  the  component. 

1.  Silicone  Oil 

The  “Series  200”  range  of  silicone  oils  from  Dow  Coming  was  chosen  to  be  the 
most  appropriate  for  the  purposes  of  this  experimental  study.  They  were  commercially 
able  to  obtain,  that  provided  a  range  of  medium  to  high  kinematic  viscosities  (50  - 
100,000  cSt)  without  any  appreciable  change  in  surface  tension  and  density.  They  were 
easy  and  safe  to  work  with,  and  they  had  advantageous  chemical  properties.  Detailed 
information  about  these  liquids  is  presented  in  Appendix-A. 

2.  Positioning  System 

In  order  for  surface  tension  to  dominate  viscous  forces,  the  Capillary  number 
(Ca)  must  be  small.  For  medium  to  high  viscosity  silicone  oils,  a  small  Capillary 
number  (10'^  -  10'®)  requires  velocities  to  be  small  as  well.  Figure  12  presents  velocities 
and  Capillary  numbers  that  are  associated  with  different  viscosity  fluids. 

The  positioning  system  control  unit  was  chosen  to  provide  these  small  velocities. 
The  detailed  information  about  the  positioning  system  is  presented  in  Appendix-B. 
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Figure  12.  Capillary  number  vs.  Velocity  plot. 


3.  Load  CeU 

From  previous  work  conducted  by  (Isik,  1999),  the  range  of  forces  was  already 
known.  Concerns  that  were  critical  in  determining  the  type  of  load  cell  were,  sensitivity, 
minimum  threshold  level,  resolution  level,  possible  drifting  in  measurements  with  time, 
signal  ou^ut,  size,  accuracy,  and  ability  of  being  used  in  tension.  The  specifications  and 
dimensions  of  the  load  cells  are  presented  in  Appendix-C. 
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4.  Optic  Systems 

Optical  system  of  the  experimental  set-up  is  composed  of  a  microscope  that  has  a 
lens  magnification  of  SOX  and  one  CCD  camera  coimected  to  a  monitor  and/or  a  VCR. 
The  CCD  camera  was  positioned  to  achieve  a  bottom  view  of  the  interfacial  liquid  layer. 

C.  EXPERIMENTAL  PROCEDURE 

Experiments  have  been  conducted  for  50,  100,  200,  10,000,  60,000,  and  100,000 
cSt  silicone  oils.  For  each  silicone  oil  viscosity  a  large  range  of  values  were  used  for  the 
initial  liquid  layer  thickness  and  the  separation  velocity.  A  brief  list  of  the  variables  and 
ranges  are  listed  below. 

LIQUID  :  Silicone  oil,  constant  surface  tension  &  density 

LIQUID  VISCOSITY  (u)  :  50, 1 00, 200, 1 0,000, 60,000, 1 00,000  cSt 

LIQUID  LAYER  THICKNESS  (L)  :  25, 50, 100  (125),  200  (250),  312.5,  500  pm 
SEPARATION  VELOCITY  (V)  :  1 .2  pm/s  -  45.5  pm/s 

Each  experiment  was  conducted  according  to  the  following  procedure  and 
schematic  diagram  of  the  apparatus  is  provide  in  Figure  13: 

1 )  Experimental  set-up  was  turned  on  to  warm-up. 

2)  Computer  program  was  started  and  control  was  transferred  to  the 
keyboard. 

3)  The  plexi-glass  disks  were  separated  to  introduce  silicone  oil. 

4)  The  silicone  oil  was  spread  on  the  lower  disk  by  an  injector. 
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5)  The  actuator  was  moved  to  position  the  plexi-glass  disks  to  the  desired 
initial  liquid  layer  thickness.  It  was  made  sure  that  there  was  no  air  bubbles  in  the 
liquid,  and  that  the  liquid  had  spread  completely  on  both  disk  surfaces. 

6)  The  velocity,  step  size  and  iterations  were  inputted  into  the  computer 
program.  The  computer  recorded  and  displayed  step  number,  time,  displacement 
and  force  in  volts  as  output.  Visual  data  could  be  monitored  or  recorded  through 
the  CCD  camera. 

7)  Each  run  was  carried  through  until  separation  of  the  liquid  layer  occurred. 

8)  These  data  have  been  analyzed  by  using  plots  of  the  force-displacement 
data  force-Capillary  number  and  the  visual  recorded  data 


-  Stationary  Top  Support 

Load  Cell 

Microscope  (SOX) 

—  Plexi-Glass  Disks 

Inward  Radial  Velocity  (Vf) 
-  Movable  Bottom  Support 


CCD  Camera 


Downward  Velocity  (V) 


Figure  13.  Schematic  Diagram  of  the  Experimental  Set-up. 
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V.  RESULTS  AND  DISCUSSION 


A.  INTRODUCTION 

The  experiments  were  carried  out  for  a  wide  range  of  viscosities.  A  wide  range  of 
separation  velocities  and  initial  liquid  layer  thicknesses  were  also  tested.  Viscosities 
ranged  from  50  cSt  to  100,000  cSt,  separation  velocities  from  1  to  47.5  pm/s,  and  initial 
liquid  layer  thicknesses  from  25  to  500  pm  were  used.  Appendix  D  presents  the 
combinations  of  these  variables. 

B.  FORCE-DISPLACEMENT  PLOTS 

The  following  plots  yield  the  resultant  force  as  the  liquid  layer  is  being  separated 
at  a  constant  velocity.  The  force  increases  almost  linearly  to  a  peak  value  (Fmax)  and 
occurs  at  a  plate  separation  distance  (dmax)  in  the  liquid  layer  as  explained  earlier.  The 
force  then  drops  rapidly  and  almost  instantly  to  a  near  zero  (steady-state)  value  and 
remains  there  with  no  further  appreciable  change.  Visually,  the  sudden  decrease  in  force 
coincides  in  the  onset  of  viscous  fingering.  The  medium  viscosity  fluids  (50  -  200  cSt) 
did  not  show  viscous  fingering  to  the  same  extent  of  the  high  viscosity  fluids  (10,000  - 
100,000),  so  the  decline  after  the  peak  value  was  more  subtle. 

After  reviewing  the  force-displacement  plots,  Fmax  and  dmax  were  dependent  on 
the  initial  liquid  layer  thickness  (lo),  the  separation  velocity  (V),  the  liquid  viscosity  (u), 
and  of  course  surface  tension  (a  special  feature  of  these  silicone  oils  is  that  the  surface 
tensions  is  relatively  constant  throughout  the  viscosities).  Each  one  of  these  variables 
were  varied  in  order  to  achieve  a  relationship  for  Fmax  and  dmax-  The  following  figures 
show  the  dependence  of  each  of  the  variables  for  both  the  medium  and  high  viscosity 
silicone  oils. 
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force  (Ibf) 


1)  Fmax,  dmax  ocu 

For  a  given  initial  liquid  layer  thickness  (lo)  and  separation  velocity  (V),  peak 
force  (Fmax)  and  the  location  of  the  peak  force  (dmax),  increase  with  increasing  liquid 
viscosity  (o)  as  shown  in  Figures  14  and  15.  The  magnitude  of  the  force  differs  greatly 
between  the  high  and  medium  viscosity  silicone  oils.  This  result  may  be  interpreted  as 
the  extent  of  the  onset  of  viscous  fingering  in  the  fluid  yields  greater  forces. 


displacement  (mm) 

Figure  14.  Force-Displacement  plot  that  shows  the  effect  of  medium  viscosity  fluids  on  and  d^ax  The 
Initial  Liquid  Layer  Thickness  (U)  is  50pm,  and  Separation  Velocity  (V)  is  24pm/s. 
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Figure  15.  Force-Displacement  plot  that  shows  the  effect  of  high  viscosity  fluids  on  Fmax  and  d^ax  The 
Initial  Liquid  Layer  Thickness  (lo)  is  312|im,  and  Separation  Velocity  (V)  is  12pni/s. 


2)  Fmax,  dmax  “cV 

For  a  given  liquid  viscosity  (u)  and  initial  liquid  layer  thickness  (lo),  peak  force 
(Fmax)  and  the  location  of  the  peak  force  (dmax),  increase  with  increasing  separation 
velocity  (Figures  16  and  17). 
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displacement  (mm) 

Figure  17.  Force-Displacement  plot  that  shows  the  effect  of  Separation  Velocity  on  F^ax  and  dm,*  The 
Initial  Liquid  Layer  Thickness  (lo)  is  500pm  and  Viscosity  (o)  is  100,000  cSt. 
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3)  Fmax,  d^ax  “c  l/lo 


For  a  given  liquid  viscosity  (u)  and  separation  velocity  (V),  peak  force  (Fmax)  and 
the  location  of  the  peak  force  (dmax),  increase  with  decreasing  initial  liquid  layer  thickness 
(lo)  (Figures  18  and  19). 
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Figure  18.  Force-Displacement  plot  that  shows  the  effect  of  Initial  Liquid  Layer  Thickness  on  and  d„ax 
Separation  Velocity  (V)  is  24pni/s  and  Viscosity  (u)  is  50  cSt. 
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4)  A  Physical  Model 

A  model  of  the  force-displacement  behavior  in  the  separation  mechanics  of  the 
layer  can  be  developed  based  on  a  creeping  flow  analysis  of  the  governing  Navier-Stokes 
equations  in  axisymmetric  cylindrical  coordinates.  The  slow  motion  allows  a  quasi¬ 
steady  analysis  in  which  viscous  effects  are  balanced  by  pressure  gradients  while  inertia 
effects  are  negligible.  The  analysis  follows  the  model  problem  suggested  by  Batchelor 
(1967)  or  Landau  &  Lifshitz  (1987)  and  indicates  that  the  separating  motion  of  the  two 
disks  creates  a  negative  suction  pressure  in  the  liquid  layer  gap  thus  resulting  a  resisting 
clamping  effect  opposite  to  the  direction  of  separation.  This  resistive  force  is  proportional 
to  liquid  viscosity  and  separation  velocity,  and  inversely  proportional  to  the  cube  of  the 
liquid  layer  thickness  at  any  instant,  a  feature  that  is  qualitatively  borne  out  by  the  current 
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study.  However  a  more  thorough  corroboration  between  model  and  experimental  data 
will  be  considered  in  future  work  once  more  extensive  data  has  been  gathered. 

C.  VISCOUS  FINGERING 

It  was  proven  that  at  the  instant  viscous  fingering  occurs,  the  interface  in  no 
longer  stable.  The  interface  motion  is  driven  by  the  separation  of  the  two  disks.  Surface 
tension  serves  as  a  stabilizing  entity  between  the  two  immiscible  fluids.  From  the 
discussion  in  the  early  chapter  on  the  Saffinan-Taylor  instability,  a  critical  wavelength 
was  determined  at  the  onset  of  viscous  fingering.  At  the  onset  of  viscous  fingering  a 
number  of  disturbances  develop  along  the  interface  called  corrugations.  At  this  point 
the  surface  tension  can  no  longer  stabilize  the  interface  and  it  is  no  longer  undistorted  and 
circular.  The  interface  is  now  unstable  and  continues  to  grow  tmstable.  If  the  instability 
proceeds  the  interface  is  no  longer  balancing  the  difference  in  viscosities  and  pressure 
differentials  to  surface  tension,  the  number  of  corrugations  has  reached  its  maximum 
value  and  coincides  with  the  maximum  resultant  force.  From  the  Saffinan-Taylor 
instability  analysis  of  the  experimental  radial  Hele-Shaw  the  maximum  wavelength  and 
number  of  corrugations  were  derived.  The  maximum  wavelength  is  the  largest 
wavelength  that  gives  rise  to  the  viscous  fingers  that  grow  the  fastest.  In  order  to  prevent 
viscous  fingering  from  occurring  and  leading  to  an  unstable  interface  the  critical 
wavelength  critical  wavelength  becomes  an  important  measuring  device  of  instability. 

A  comparison  of  various  separation  velocities  and  initial  liquid  layer  thickness  are 
considered  below.  The  maximum  wavelength  and  number  of  corrugations  for  each  run 
was  calculated  and  compared  to  the  picture  at  the  same  instant.  Table  1  contains  the 
calculated  quantities. 
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Kinematic 
Viscosity  (cSt) 

10000 

10000 

10000 

10000 

0.021 

0.021 

0.021 

0.021 

Radius  (m) 

0.0254 

Initial  Liquid  Layer 
Thickness  (m) 

0.0002 

0.0003125 

0.0003125 

Displacement  (mm) 

0.029 

0.051 

0.0328 

0.051 

Radial  Velocity 
(m/s) 

3.75x10-^ 

2.24x10'^ 

3.0x10-^ 

Separation  Velocity 
(m/s) 

24x10*^ 

12x10-'’ 

24x10-'’ 

Maximum 
Wavelength  (m) 

ma 

Number  of 
Corrugations 

80 

84 

47 

52 

Table  1.  Experiments  conducted  for  10,000  cSt  Silicone  Oil,  Ca! 

'Culating  Maximum  Wavelength  and 

Number  of  Corrugations. 


As  initial  liquid  layer  thickness  (lo)  decreases,  viscous  fingering  becomes  more 
pronounced,  leading  to  an  increase  in  number  of  corrugations  and  a  decrease  in  mavirmim 
wavelength.  Similarly,  separation  velocity  (V)  increases,  viscous  fingering  becomes 
more  pronounced  resulting  in  an  increase  in  the  number  of  corrugations  and  a  decrease  in 
maximum  wavelength.  Also,  separation  velocity  (V)  and  the  liquid  viscosity  (u)  also 
affect  the  degree  of  viscous  fingering.  Higher  the  separation  velocity  and  liquid 
viscosity,  greater  is  the  extent  of  viscous  fingering. 

The  following  pictures  were  taken  at  the  maximum  growth  rate,  although  the 
actual  number  of  corrugations  can  be  difficult  to  identify  the  same  trend  as  shown  above 
for  the  maximum  wavelength  and  number  of  corrugations  can  be  seen. 
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Figure  20. 10,000  cSt  Silicone  Oil,  Initial  Liquid  Layer  Thickness  of 200Mm  and  Separation  Velocity  of 
12pm/s  (Scale  of  ldiv=lmm),  N»80. 


Figure  21.  10,000  cSt  Silicone  Oil,  Initial  Liquid  Layer  Thickness  of  200nm  and  Separation  Velocity  of 
24pm/s  (Scale  of  ldiv=lniin),  N«84. 
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Figure  22.  10,000  cSt  Silicone  Oil,  Initial  Liquid  Layer  Thickness  of  3 12.5|im  and  Separation  Velocity  of 
12pm/s  (Scale  of  ldiv=lmm),  Nw47. 


Figure  23,  10,000  cSt  Silicone  Oil,  Initial  Liquid  Layer  Thickness  of  312.5pm  and  Separation  Velocity  of 
24pm/s  (Scale  of  ldiv=lnim),  N«52. 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Well-defined  trends  in  the  force-displacement  behavior  of  thin  liquid  layers  were 
identified  from  the  analysis  of  the  recorded  quantitative  and  qualitative  data  in  this 
experimental  study.  A  wide  range  of  viscosities  were  examined  and  all  force- 
displacement  curves  have  lead  to  the  observation  of  the  presence  of  a  peak  force  (Fmax) 
whose  location  and  value  depend  on  initial  liquid  layer  thickness,  separation  velocity, 
liquid  viscosity,  and  surface  tension. 

The  presence  of  the  peak  force  is  more  pronounced  and  its  occurrence  coincides 
with  the  onset  of  the  viscous  fingering  mechanism.  The  approach  to  the  peak  force  (Fmax) 
on  the  force-displacement  plot,  is  accompanied  by  an  instantaneous  and  sudden  decrease 
in  force  which  is  followed  by  a  near  zero  (steady-state)  value  with  no  further  appreciable 
change.  The  linear  stability  analysis  allows  a  fairly  accurate  prediction  of  the  onset  of 
the  viscous  fingering  behavior  and  correlates  well  with  the  experimental  data. 

Further  research  on  this  subject  is  highly  desirable  to  widen  the  parametric  study 
to  other  viscosity  ranges  and  surface  tension.  A  more  accurate  characterization  of  the 
interface  velocity  is  crucial  to  further  the  understanding  of  the  onset  of  viscous  fingering 
and  its  dependence  on  the  Capillary  number.  A  review  and  refinement  of  the  analytical 
model  may  also  be  conducted  and  balanced  with  the  necessary  boundary  conditions, 
which  could  lead  to  corroboration  with  the  different  observed  regimes. 

In  addition,  a  more  advanced  and  integrated  testing  apparatus  with  better 
measurement  abilities  of  the  interface  velocity,  layer  thickness,  actuator  motion  and 
interface  flatness  and  parallelism  is  required.  It  would  also  be  desirable  to  increase  the 
data  acquisition  rate  to  at  least  the  milli-second  level  in  order  obtain  more  accurate  force- 
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displacement  information.  The  use  of  an  optical  grade  glass,  instead  of  plexi-glass,  for 
the  disk  would  be  advantageous  to  ensure  smooth,  flat  and  parallel,  and  hard  (scratch- 
resistant)  surfaces. 
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APPENDIX -A 


SILICONE  OIL 

Medium  Viscosity  Silicone  Fluids: 

200®  Fluid,  50  cSt 
200®  Fluid,  100  cSt 
200®  Fluid,  200  cSt 

High  Viscosity  Silicone  Fluids: 

200®  Fluid,  10,000  cSt 
200®  Fluid,  60,000  cSt 
200®  Fluid,  100,000  cSt 

Description: 

200®  Fluids  from  Dow  Coming,  50  —  1,000  centistokes  (cSt)  are  medium,  and 
10,000  —  100,000  centistokes  (cSt),  are  high  viscosity  polydimethylsiloxane  polymers 
manufactured  to  yield  linear  polymers  with  average  kinematic  viscosities  ranging  from 
50  -  1000  cSt  and  10,000  -  100,000  cSt. 

Composition: 

Linear  polydimethylsiloxane  polymers  characteristically  have  the  following 
typical  chemical  composition: 

(CHs)sSiO[SiO(CHs)J,^i(CH0, 

Commercial  bulk  polymerized  dimethl  silicone  fluids,  such  as  200®  Fluids  typically 
contain  trace  amounts  of  process  impiirities. 

Benefits: 

200®  Fluids,  50  —  100,000  cSt,  have  the  following  product  characteristics; 

►  Clear  ►  Nongreasy  ►  Nonocclusive 

►  Non-stinging  on  skin 
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200®  Fluids,  50  -  100,000  cSt,  when  compared  with  other  materials  that  may  be 
substituted  in  a  given  application,  may  offer  one  or  more  of  these  comparative 
characteristics: 

►  High  compressibility  ►  High  damping  action 

►  High  oxidation  resistance  ►  High  shear-ability  without  breakdown 

►  High  temperature  serviceability  ►  High  compatibility 

►  High  spreadibility  High  water  repellency 

►  Low  fire  hazard  ►  Low  odor 

►  Low  reactivity  ►  Low  surface  energy 

►  Low  temperature  serviceability  ►  Low  vapor  pressure 

►  Good  heat  stability  ►  Soft  feel  and  lubricity  on  skin 

►  Good  leveling  and  easy  rubout 
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APPENDIX -B 


POSITIONING  SYSTEM 


Model  855  Programmable  Controller  System  (Newport  Corporation) 

The  Model  855  Programmable  Controller  is  a  microprocessor-based  system  that 
allows  simultaneous  direct  or  programmable  control  of  up  to  four  Newport  linear 
actuators  or  rotary  stages.  It  is  simple  but  powerful  vocabulary  of  mnemonic  commands 
allows  straightforward  programming  and  control  via  the  handheld  855K  Keypad/Display. 
The  855C  Controller: 


The  855C  Programmable  Controller  is  the  nucleus  of  a  system  that  automatically 
controls  up  to  foiir  Newport  precision  positioners  and  stages.  Its  large,  easy-to-leam 
instruction  set  and  standard  RS-232C  and  IEEE-488  interface  ports  allow  it  to  work 
closely  with  external  computers  and  other  data  devices.  Its  programmability  provides 
stand-alone  automatic  control  of  actuator  motion,  yet  no  knowledge  of  programming 
techniques  is  required.  It  also  supports  the  optional  855K  handheld  Keypad/Display  for 
convenient  data  entry,  control  and  program  editing. 


Specifications: 


Actuator  Control: 


CPU: 

RAM: 

PROM: 

EPROM: 

I/O: 

Power: 

Dimensions: 


Provides  ±15  VDC,  1.5  Amp  power  and  ±10V  velocity  control 
signals  for  up  to  4  Newport  linear  positioning  devices.  Receives 
and  decodes  dual  output,  90°  phase,  ±12.5,  -0.5  V  encoder  pulses 
and  limit  signals. 

8  bit  6809  high-performance  micro-processor 
4  Kbytes 
24  Kbytes 
2  Kbytes 

RS-232C  Serial  Port/IEEE-488  Parallel  Port 

110-220  VAC,  50/60  Hz  and  850  nxA  with  all  actuators  Mly 

loaded. 

5.25”  high  X  17”  wide  X  12”  deep 
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The  855K  Keyboard/Display: 

The  855K  can  be  used  to  completely  control  and  program  the  855C.  It  is  handy 
for  manual  control  of  actuator  position  and  for  the  reviewing  and  editing  855C  programs. 
Its  backlit  alphanumeric  liquid-crystal  display  shows  position  information  for  the  four 
drivers,  programming  instructions  for  the  855C,  and  messages  from  the  855C.  The 
855K’s  keypad  has  35  durable,  tactile-response  membrane  keys  for  command  and  data 
entry.  The  855K  connects  to  the  855C  via  a  coiled  cable. 

The  Linear  Actuators: 

The  855C  is  compatible  with  Newport’s  850  Series  Linear  Actuators.  Each  of 
these  devices  use  DC  motor  drives  with  integral  optical  encoders  for  smooth  operation 
and  high  resolution.  The  855C  fully  supports  their  resolution  and  range.  Only  one 
actuator  is  used  for  this  experimental  study. 
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APPENDIX -C 


LOAD  CELL 

Since  the  viscosities  of  the  fluids  differed  greatly  from  50  cSt  to  100,000  cSt, 
there  was  a  need  to  use  two  different  load  cells  for  the  experiments.  For  higher 
viscosities  an  MLP-100  (100  lb  maximum  load)  and  for  the  lower  viscosities  the  MDP- 
2.5  (2.5  lb  maximum  load)  were  used.  The  specifications  are  listed  below  for  each  load 
cell. 


MLP-100  Specifications: 

MmimpRQFiiEmm 

mmmL  /  fEmmm  mispmmH 

The  model  MLP  Series  load  cells  were  designed  with  economy  as  first  priority. 
The  MLP  Series  are  anodized  aluminum  with  a  unique  low  profile  design,  which 
provides  excellent  stability  for  in  line  application  for  tension  and/or  compression. 


f  la  4  COimiaOR  COLOR  COOED. 

I  SHiElDEPCAeLcSUffUEp 

tHREAD  TYPICAL 
SOTHE^DS 


43 


^eimmoHS 


Rated  Ou^xit  {R.O.}; 
Noiilinearit)^ 
Mfyi^resis: 
Nonrep^tafailHy: 

Comper^!todiWm{>.  Range: 

1^p.  Rteige: 
Effi^dn  Outp^ 
Terdp.  0dct  Zero: 
Resistenoe: 
^itettonV^ttage: 

4'- ..VWsigM: 


2  mV/V  nortiinid 
O.O^df  RlO. 
0.0S%ofR-O. 
O.O^idfR.O. 
1.0%  Of  R  04 
60*  to  160# 
^'to^OT 
0.005%  ofiLoad^ 
0i>05%ofR;O7# 
360  phrns  noMinal 
10VDC 
150%ofR.O: 
1qz.«R  ranges 


DIMENSIONS  (INCHES) 


CAPACnV  1>USA0 

MODEL  LBS.  L  LI  W  W1  H  THREAD  DEPTH 


lilLP-IO 

10 

14604 

'-^4eoo."::^:’^:3754-  js 

4vl0^4{>: 

.;-JS00 

25 

HOMO 

50 

'  12304^4:; 

MLP>7S 

75 

ltLP*10Q 

100 

1J664 

msEm 

iaF-150 

ISO 

ML^200 

200 

4  ^  ^  J50 

3«^4 

'jsrs- .A 

jyojxsoo 

900 

500 

t4;;sfio^-'^^;i:v:55o "  -1:00 ' 

,37s 

IILP-750 

750 

IIUMK 

1,000 

^  7^  1,00 

3^ 

E» 
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MDB-2.5  Specifications: 


The  MDB  Series  was  designed  to  help  fill  the  growing  need  for  a  greater  selection 
of  high  accuracy  load  cells  for  use  in  space  limited  applications.  The  anodized  aluminum 
MDB's  are  compliant  tension  and  compression,  therefore,  a  good  choice  for  in  line 
through  zero  applications,  as  well  as  single  direction  tension  or  compression. 


\.U4-2aUNF 
,2SO  DEEP  BOTH  ENt® 


SPESlFiCATlOMS 
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APPENDIX -D 


EXPERIMENTAL  DATA 


50, 100^  200 
cSt 

Silicone  Oil 


Table  1.  50, 100, 200  cSt  silicone  oi 


Initial  Liquid  Layer 
Thickness  (L) 

25  pm 

48 

24 

Velocity  12 

( LLin/s^ 

VM'  J  2  4 

1.2 

Initial  Liquid  Layer 
Thickness  (L) 
50pm 

48 

24 

Velocity  12 

(ijm/s)  4.5 

^  2.4 

1.2 

Initial  Liquid  Layer 
Thickness  (L) 
100pm 

48 

24 

Velocity  _  12 

('Lim/s'l 

^  2  4 

1.2 

Initial  Liquid  Layer 
Thickness  (L) 

200pm 

48 

24 

Velocity  12 

(|4m/s)  4.5 

2.4 

1.2 

Initial  Liquid  Layer 
Thickness  (L) 
312pm 

48 

24 

Velocity  12 

lum/s)  4.5 

2.4 

1.2 

Initial  Liquid  Layer 
Thickness  (L) 
500pm 

48 

24 

Velocity  12 

(um/s)  4.5 

2.4 

1.2 

s  and  their  varying  parameters. 
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Initial  Liquid  Layer 
Thickness  (L) 
125|mi 

48 

24 

Velocity  12 

(Min/s)  - II - 

1.2 

10,000, 

60,000, 

Initial  Liquid  Layer 
Thickness  (I,) 
250|jm 

48 

24 

Velocity  12 

(jLim/s) 

2.4 

1.2 

100,000  cSt 
Silicone  Oil 

Initial  Liquid  Layer 
Thickness  (L) 
312pm 

48 

24 

Velocity  12 

(pm/s)  - ^ - 

1.2 

Initial  Liquid  Layer 
Thickness  (lo) 
500pm 

48 

24 

Velocity  j2 

(pm/s)  4.5 

2.4 

1.2 

Table  1.  10,000,  60,000, 100,000  cSt  silicone  oils  and  their  varying  parameters. 
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